Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 3D computation of incompressible multi-phase flows described by conservative phase field models. We focus here on the case of density matched fluids with different viscosity (Model H). The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow’s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow induced drop coalescence.